\(\int \frac {e^{-a-b x} (a+b x)^3}{x} \, dx\) [60]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [A] (verified)
   Fricas [A] (verification not implemented)
   Sympy [A] (verification not implemented)
   Maxima [A] (verification not implemented)
   Giac [A] (verification not implemented)
   Mupad [B] (verification not implemented)

Optimal result

Integrand size = 21, antiderivative size = 102 \[ \int \frac {e^{-a-b x} (a+b x)^3}{x} \, dx=-2 e^{-a-b x}-3 a e^{-a-b x}-3 a^2 e^{-a-b x}-2 b e^{-a-b x} x-3 a b e^{-a-b x} x-b^2 e^{-a-b x} x^2+a^3 e^{-a} \operatorname {ExpIntegralEi}(-b x) \]

[Out]

-2*exp(-b*x-a)-3*a*exp(-b*x-a)-3*a^2*exp(-b*x-a)-2*b*exp(-b*x-a)*x-3*a*b*exp(-b*x-a)*x-b^2*exp(-b*x-a)*x^2+a^3
*Ei(-b*x)/exp(a)

Rubi [A] (verified)

Time = 0.10 (sec) , antiderivative size = 102, normalized size of antiderivative = 1.00, number of steps used = 9, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.190, Rules used = {2230, 2225, 2209, 2207} \[ \int \frac {e^{-a-b x} (a+b x)^3}{x} \, dx=e^{-a} a^3 \operatorname {ExpIntegralEi}(-b x)-3 a^2 e^{-a-b x}-b^2 x^2 e^{-a-b x}-3 a e^{-a-b x}-3 a b x e^{-a-b x}-2 e^{-a-b x}-2 b x e^{-a-b x} \]

[In]

Int[(E^(-a - b*x)*(a + b*x)^3)/x,x]

[Out]

-2*E^(-a - b*x) - 3*a*E^(-a - b*x) - 3*a^2*E^(-a - b*x) - 2*b*E^(-a - b*x)*x - 3*a*b*E^(-a - b*x)*x - b^2*E^(-
a - b*x)*x^2 + (a^3*ExpIntegralEi[-(b*x)])/E^a

Rule 2207

Int[((b_.)*(F_)^((g_.)*((e_.) + (f_.)*(x_))))^(n_.)*((c_.) + (d_.)*(x_))^(m_.), x_Symbol] :> Simp[(c + d*x)^m*
((b*F^(g*(e + f*x)))^n/(f*g*n*Log[F])), x] - Dist[d*(m/(f*g*n*Log[F])), Int[(c + d*x)^(m - 1)*(b*F^(g*(e + f*x
)))^n, x], x] /; FreeQ[{F, b, c, d, e, f, g, n}, x] && GtQ[m, 0] && IntegerQ[2*m] &&  !TrueQ[$UseGamma]

Rule 2209

Int[(F_)^((g_.)*((e_.) + (f_.)*(x_)))/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[(F^(g*(e - c*(f/d)))/d)*ExpInteg
ralEi[f*g*(c + d*x)*(Log[F]/d)], x] /; FreeQ[{F, c, d, e, f, g}, x] &&  !TrueQ[$UseGamma]

Rule 2225

Int[((F_)^((c_.)*((a_.) + (b_.)*(x_))))^(n_.), x_Symbol] :> Simp[(F^(c*(a + b*x)))^n/(b*c*n*Log[F]), x] /; Fre
eQ[{F, a, b, c, n}, x]

Rule 2230

Int[(F_)^((c_.)*(v_))*(u_)^(m_.)*(w_), x_Symbol] :> Int[ExpandIntegrand[F^(c*ExpandToSum[v, x]), w*NormalizePo
werOfLinear[u, x]^m, x], x] /; FreeQ[{F, c}, x] && PolynomialQ[w, x] && LinearQ[v, x] && PowerOfLinearQ[u, x]
&& IntegerQ[m] &&  !TrueQ[$UseGamma]

Rubi steps \begin{align*} \text {integral}& = \int \left (3 a^2 b e^{-a-b x}+\frac {a^3 e^{-a-b x}}{x}+3 a b^2 e^{-a-b x} x+b^3 e^{-a-b x} x^2\right ) \, dx \\ & = a^3 \int \frac {e^{-a-b x}}{x} \, dx+\left (3 a^2 b\right ) \int e^{-a-b x} \, dx+\left (3 a b^2\right ) \int e^{-a-b x} x \, dx+b^3 \int e^{-a-b x} x^2 \, dx \\ & = -3 a^2 e^{-a-b x}-3 a b e^{-a-b x} x-b^2 e^{-a-b x} x^2+a^3 e^{-a} \text {Ei}(-b x)+(3 a b) \int e^{-a-b x} \, dx+\left (2 b^2\right ) \int e^{-a-b x} x \, dx \\ & = -3 a e^{-a-b x}-3 a^2 e^{-a-b x}-2 b e^{-a-b x} x-3 a b e^{-a-b x} x-b^2 e^{-a-b x} x^2+a^3 e^{-a} \text {Ei}(-b x)+(2 b) \int e^{-a-b x} \, dx \\ & = -2 e^{-a-b x}-3 a e^{-a-b x}-3 a^2 e^{-a-b x}-2 b e^{-a-b x} x-3 a b e^{-a-b x} x-b^2 e^{-a-b x} x^2+a^3 e^{-a} \text {Ei}(-b x) \\ \end{align*}

Mathematica [A] (verified)

Time = 0.18 (sec) , antiderivative size = 52, normalized size of antiderivative = 0.51 \[ \int \frac {e^{-a-b x} (a+b x)^3}{x} \, dx=e^{-a-b x} \left (-2-3 a^2-2 b x-b^2 x^2-3 a (1+b x)+a^3 e^{b x} \operatorname {ExpIntegralEi}(-b x)\right ) \]

[In]

Integrate[(E^(-a - b*x)*(a + b*x)^3)/x,x]

[Out]

E^(-a - b*x)*(-2 - 3*a^2 - 2*b*x - b^2*x^2 - 3*a*(1 + b*x) + a^3*E^(b*x)*ExpIntegralEi[-(b*x)])

Maple [A] (verified)

Time = 0.21 (sec) , antiderivative size = 96, normalized size of antiderivative = 0.94

method result size
meijerg \({\mathrm e}^{-a} \left (2-\frac {\left (3 b^{2} x^{2}+6 b x +6\right ) {\mathrm e}^{-b x}}{3}\right )+3 \,{\mathrm e}^{-a} a \left (1-\frac {\left (2 b x +2\right ) {\mathrm e}^{-b x}}{2}\right )+3 \,{\mathrm e}^{-a} a^{2} \left (1-{\mathrm e}^{-b x}\right )+{\mathrm e}^{-a} a^{3} \left (\ln \left (x \right )+\ln \left (b \right )-\ln \left (b x \right )-\operatorname {Ei}_{1}\left (b x \right )\right )\) \(96\)
risch \(-b^{2} {\mathrm e}^{-b x -a} x^{2}-a^{3} {\mathrm e}^{-a} \operatorname {Ei}_{1}\left (b x \right )-3 a b \,{\mathrm e}^{-b x -a} x -3 a^{2} {\mathrm e}^{-b x -a}-2 b \,{\mathrm e}^{-b x -a} x -3 a \,{\mathrm e}^{-b x -a}-2 \,{\mathrm e}^{-b x -a}\) \(97\)
derivativedivides \(-a^{2} {\mathrm e}^{-b x -a}+a \left (\left (-b x -a \right ) {\mathrm e}^{-b x -a}-{\mathrm e}^{-b x -a}\right )-\left (-b x -a \right )^{2} {\mathrm e}^{-b x -a}+2 \left (-b x -a \right ) {\mathrm e}^{-b x -a}-2 \,{\mathrm e}^{-b x -a}-a^{3} {\mathrm e}^{-a} \operatorname {Ei}_{1}\left (b x \right )\) \(113\)
default \(-a^{2} {\mathrm e}^{-b x -a}+a \left (\left (-b x -a \right ) {\mathrm e}^{-b x -a}-{\mathrm e}^{-b x -a}\right )-\left (-b x -a \right )^{2} {\mathrm e}^{-b x -a}+2 \left (-b x -a \right ) {\mathrm e}^{-b x -a}-2 \,{\mathrm e}^{-b x -a}-a^{3} {\mathrm e}^{-a} \operatorname {Ei}_{1}\left (b x \right )\) \(113\)

[In]

int(exp(-b*x-a)*(b*x+a)^3/x,x,method=_RETURNVERBOSE)

[Out]

exp(-a)*(2-1/3*(3*b^2*x^2+6*b*x+6)*exp(-b*x))+3*exp(-a)*a*(1-1/2*(2*b*x+2)*exp(-b*x))+3*exp(-a)*a^2*(1-exp(-b*
x))+exp(-a)*a^3*(ln(x)+ln(b)-ln(b*x)-Ei(1,b*x))

Fricas [A] (verification not implemented)

none

Time = 0.27 (sec) , antiderivative size = 50, normalized size of antiderivative = 0.49 \[ \int \frac {e^{-a-b x} (a+b x)^3}{x} \, dx=a^{3} {\rm Ei}\left (-b x\right ) e^{\left (-a\right )} - {\left (b^{2} x^{2} + {\left (3 \, a + 2\right )} b x + 3 \, a^{2} + 3 \, a + 2\right )} e^{\left (-b x - a\right )} \]

[In]

integrate(exp(-b*x-a)*(b*x+a)^3/x,x, algorithm="fricas")

[Out]

a^3*Ei(-b*x)*e^(-a) - (b^2*x^2 + (3*a + 2)*b*x + 3*a^2 + 3*a + 2)*e^(-b*x - a)

Sympy [A] (verification not implemented)

Time = 3.43 (sec) , antiderivative size = 70, normalized size of antiderivative = 0.69 \[ \int \frac {e^{-a-b x} (a+b x)^3}{x} \, dx=\left (a^{3} \operatorname {Ei}{\left (- b x \right )} - 3 a^{2} e^{- b x} - 3 a \left (b x e^{- b x} + e^{- b x}\right ) - b^{2} x^{2} e^{- b x} - 2 b x e^{- b x} - 2 e^{- b x}\right ) e^{- a} \]

[In]

integrate(exp(-b*x-a)*(b*x+a)**3/x,x)

[Out]

(a**3*Ei(-b*x) - 3*a**2*exp(-b*x) - 3*a*(b*x*exp(-b*x) + exp(-b*x)) - b**2*x**2*exp(-b*x) - 2*b*x*exp(-b*x) -
2*exp(-b*x))*exp(-a)

Maxima [A] (verification not implemented)

none

Time = 0.23 (sec) , antiderivative size = 69, normalized size of antiderivative = 0.68 \[ \int \frac {e^{-a-b x} (a+b x)^3}{x} \, dx=a^{3} {\rm Ei}\left (-b x\right ) e^{\left (-a\right )} - 3 \, {\left (b x + 1\right )} a e^{\left (-b x - a\right )} - 3 \, a^{2} e^{\left (-b x - a\right )} - {\left (b^{2} x^{2} + 2 \, b x + 2\right )} e^{\left (-b x - a\right )} \]

[In]

integrate(exp(-b*x-a)*(b*x+a)^3/x,x, algorithm="maxima")

[Out]

a^3*Ei(-b*x)*e^(-a) - 3*(b*x + 1)*a*e^(-b*x - a) - 3*a^2*e^(-b*x - a) - (b^2*x^2 + 2*b*x + 2)*e^(-b*x - a)

Giac [A] (verification not implemented)

none

Time = 0.29 (sec) , antiderivative size = 95, normalized size of antiderivative = 0.93 \[ \int \frac {e^{-a-b x} (a+b x)^3}{x} \, dx=-b^{2} x^{2} e^{\left (-b x - a\right )} + a^{3} {\rm Ei}\left (-b x\right ) e^{\left (-a\right )} - 3 \, a b x e^{\left (-b x - a\right )} - 3 \, a^{2} e^{\left (-b x - a\right )} - 2 \, b x e^{\left (-b x - a\right )} - 3 \, a e^{\left (-b x - a\right )} - 2 \, e^{\left (-b x - a\right )} \]

[In]

integrate(exp(-b*x-a)*(b*x+a)^3/x,x, algorithm="giac")

[Out]

-b^2*x^2*e^(-b*x - a) + a^3*Ei(-b*x)*e^(-a) - 3*a*b*x*e^(-b*x - a) - 3*a^2*e^(-b*x - a) - 2*b*x*e^(-b*x - a) -
 3*a*e^(-b*x - a) - 2*e^(-b*x - a)

Mupad [B] (verification not implemented)

Time = 0.23 (sec) , antiderivative size = 69, normalized size of antiderivative = 0.68 \[ \int \frac {e^{-a-b x} (a+b x)^3}{x} \, dx=-{\mathrm {e}}^{-a-b\,x}\,\left (b^2\,x^2+2\,b\,x+2\right )-3\,a^2\,{\mathrm {e}}^{-a-b\,x}-3\,a\,{\mathrm {e}}^{-a-b\,x}\,\left (b\,x+1\right )-a^3\,{\mathrm {e}}^{-a}\,\mathrm {expint}\left (b\,x\right ) \]

[In]

int((exp(- a - b*x)*(a + b*x)^3)/x,x)

[Out]

- exp(- a - b*x)*(2*b*x + b^2*x^2 + 2) - 3*a^2*exp(- a - b*x) - 3*a*exp(- a - b*x)*(b*x + 1) - a^3*exp(-a)*exp
int(b*x)